Cotunneling current through quantum dots with phonon-assisted spin-flip processes 
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We consider cotunneling through a quantum dot in the presence of spin-flip processes induced by 
the coupling to acoustic phonons of the surrounding. An expression for the phonon-assisted cotun- 
neling current is derived by means of a generalized Schrieffer- Wolff transformation. The influence 
of the spin-phonon coupling on the heating of the dot is considered. The result is evaluated for the 
case of a parabolic semiconductor quantum dot with Rashba and Dresselhaus spin-orbit coupling 
and a method for the determination of the spin-phonon relaxation rate is proposed. 
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I. INTRODUCTION 

In the past years, prospective applications in spintron- 
ics and quantum computationiS*^ have lead to a grow- 
ing interest in the magnetic properties of nanoscale sys- 
tems. A variety of systems ranging from semiconduc- 
tor quantum dots over molecular magnets^SiSiS down to 
single atoms on surfaces have been studied in detailpi^^ 
For the investigation of such systems, transport measure- 
ments provide an excellent tool. Recently, the Zeeman 
splitting of two spin-levels has been measured via inelas- 
tic tunneling spectroscopy, both in GaAs quantum dots 12 
and in single Mn atoms on a surface which have been ad- 
dressed by an STM tip«i°i Using the same technique, the 
singlet-triplet splitting in few-electron quantum dots has 
been determined. 13 

For the applications mentioned at the beginning it is 
crucial that effects like relaxation and the loss of quantum 
coherence resulting from the coupling of the nanosystem 
to its surrounding environment are sufficiently weak. In 
particular, a detailed knowledge about these effects is a 
prerequisite for the assessment of the suitability of a spe- 
cific system for application purposes. In the case of the 
electron spin in semiconductor quantum dots, relaxation 
at low magnetic fields is dominated by the hyperfine cou- 
pling to the nuclei jiii^iiii while at higher fields the spin- 
orbit interaction induced coupling to acoustical phonons 
is most relevant piLiS It turns out that both effects are suf- 
ficiently weak, and correspondingly long spin-relaxation 
times have been reportedpis^SiSi For the single-atom ex- 
periments reported in Ref. IToL the strength of different 
coupling mechanisms is not yet clear, however, and the 
measurement schemes which were used in the semicon- 
ductor case are not directly applicable. 

In the present work, we study to what extent inelas- 
tic cotunneling spectroscopy can yield information about 
intrinsic spin-flip processes due to a coupling to bosonic 
degrees of the freedom in the surrounding. We put for- 
ward a general approach for the calculation of the co- 
tunneling current* 2 ^ 3 , in the presence of such spin-flip 
processes. Subsequently, we apply the general formalism 
to the case of cotunneling across GaAs/AlGaAs semi- 
conductor quantum dots and propose a method for the 



determination of the spin-relaxation rate. 

The outline of the present work is as follows. In Sec.HTl 
we start by introducing the model for the quantum dot 
coupled via a tunneling Hamiltonian to two fermionic 
leads and via a spin-phonon coupling to a bosonic envi- 
ronment. In the next section, we put forward a general- 
ized Schrieffer- Wolff transformation, which, in the limit 
of a weak spin-phonon coupling, allows one to eliminate 
to lowest order the dot-lead coupling. We then derive, in 
Sec. lIVI expressions for the elastic and inelastic cotunnel- 
ing current in the presence of the spin-phonon coupling. 
For the evaluation of the current, the knowledge of the 
occupation probabilities of the spin-states on the dot is 
required. Their dynamics is governed by a master equa- 
tion, which is discussed in Sec. The experimentally 
most relevant quantity, the differential conductance as a 
function of the bias voltage, is evaluated in Sec. IVII In 
Sec. lVIII we apply the general result to a realistic model 
describing a recent experiment by Kogan et al. (Ref. [13). 
Sec. VIII briefly discusses the low-temperature behavior. 
The final conclusions are presented in Sec. IIXI 



II. MODEL 

We consider transport across a quantum dot in a spin- 
1/2 ground-state contacted to two leads in the presence 
of an intrinsic spin-flip mechanism due to the coupling of 
the dot to a phonon bath (see Fig. The total system 
consisting of dot, leads and phonons is described by the 
Hamiltonian 



H = Hq + Hi 



(1) 



Here, Hq is the Hamiltonian of the isolated dot, the ideal 
leads I = L, R and the phonons, i.e., 

H Q = ^2 E a n a +U n T n x +y^ Q k ni kg .+y^ faj q n q . (2) 

<T=TJ ika q 

The number operator of the dot spin is given by n a — 
d\d a , where d G (rfj.) destroys (creates) an electron with 
spin a =t,| on the dot. An externally applied mag- 
netic field B z in z-direction leads to a Zeeman splitting 
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(4) 



So far, except for the presence of two instead of one lead, 
the Hamiltonian corresponds to the one of the single-site 
Anderson model^- which has been extensively studied in 
the literature^ 

Finally, we assume that the dot spin interacts with the 
phonon system via a spin-phonon coupling of the form 



(5) 



FIG. 1: Sketch of the quantum dot coupled via the tunneling 
Hamiltonian Ht to two leads at electro-chemical potentials /il 
and ^ir and to a phonon bath via the spin-phonon coupling 
Hamiltonian Hs- The dot can be either empty (0), in one of 
the spin-directions f and J., which are split by the Zeeman 
energy Aj.j, or in the singlet state S. 



In order to guarantee the hermiticity of the coupling 
Hamiltonian, the coefficients have to fulfil M q i = M^ q i 
for all q and i = x,y. Furthermore, we have intro- 
duced the spin operators u x = djd-j- + djdj. and o y = 

— d^di). Again, we can characterize the coupling 
by a spectral density, which is defined as 



Mr 



— i?j = —g\i-&B z of the single-particle dot lev- 
els. Furthermore, due to electron-electron interaction, 
the simultaneous occupation of the dot by two electrons 
costs the additional energy U. The leads I — L, R are 
modeled as Fermi liquids, i.e., by non-interacting quasi- 
particles with energy Qk, where k is the wave- vector and 
<t the spin. The number operator n^ko- = cj kcr Qko- for the 
corresponding state is defined in terms of the fermionic 
creation and destruction operators c| ko . and Qk<r, respec- 
tively. Each lead is filled up to an electro-chemical poten- 
tial fj,#, and thus the mean occupation numbers obey the 
Fermi distributions fe{e) = {exp [(e — ^e)/k^T\ + 1} . 
Here, an externally applied bias voltage V maps to a dif- 
ference — Mr = eV . In the present work, we consider 
transport across a singly-occupied dot in the cotunncling 
regime. This imposes the condition E a -C fig <C E a + U 
for cr =f,l and I = L, R. The last term in Eq. © 
describes the phonons of the environment, which con- 
sist of a collection of modes characterized by a wave- 
vector q with occupation numbers n q = a q a q . Here, 
the bosonic operators a q (a q ) create (destroy) a phonon 
in the mode q. Note that for reasons of notational 
simplicity, a possible branch index has been suppressed 
here, but can be added in the end results. At equilib- 
rium, the phonons are described by the Bose-distribution 
ra B (e) = [exp( e /fc B T) - l]- 1 . 

The coupling of the dot to the leads is modeled by the 
tunneling Hamiltonian 



H T = ^2 T ik C^da + h.c. , 



(3) 



Ikcr 



where, for reasons of simplicity, we have assumed spin- 
independent tunnel matrix elements T^. For later use, 
we introduce the corresponding tunneling rates (per spin- 



D(ui) = J ^ + iM i>v\ 2 ^ ~ w q) ■ ( 6 ) 
q 

Note that the coupling does not contain dephasing terms 
proportional to a z = dld^ — d]d± or to the total number 
of electrons n = + n± on the dot. For the electron 
spin in a quantum dot this is justified to a very good ap- 
proximation, as was shown in Rcf. In this regard, the 
present model differs, e.g., from the so-called Anderson- 
Holstein model, 26 - 2 i which only contains coupling terms 
of this kind. 



III. SCHRIEFFER- WOLFF TRANSFORMATION 

Since we are interested in transport in the cotunneling 
regime where the dot occupancy is only changed virtu- 
ally, it is advantageous to eliminate to lowest order the 
tunneling Hamiltonian Ht, which changes the number 
of electrons on the dot by one. A standard method to 
achieve this goal is the so-called Schrieffer- Wolff transfor- 
mation, 28 which we generalize to the case where a spin- 
flip coupling Hs is present. Accordingly, we perform the 
canonical transformation to the Hamiltonian 



H = e s He~ s = H 



H S + ±[S,H T ]+0(H%), (7) 



where the generator S of the transformation has to fulfil 
the condition [S, H + Hs] + Ht = 0. In order to obtain 
an explicit expression for the generator S, we now use 
that, as has been discussed in the Introduction, in all 
relevant cases, the spin-flip coupling is weak and can be 
treated perturbatively. Hence, we can expand 
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L^[L^H Tl H s ] 



0(H 



(8) 
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Here, Lq is the Liouvillian of the decoupled dot, leads and 
phonon system^ We, thus, arrive at the transformed 
Hamiltonian 

H s» H +H s + - [L 1 H T ,H T ] + -[L 1 [L Q 1 H T , Hs],Ht] . 

(9) 

The first commutator corresponds to the standard 
Schrieffer- Wolff transformation 2 ^ and yields a contribu- 
tion 

O t^O 1 ^T, Hi] = #dir,ex + Hren + H 2c ■ (10) 

Here, the first term contains the direct and exchange in- 
teraction between the dot and the leads 
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d\ d a cj kcr Cfk'ff) + h.c. (11) 



We denote flipped spins by | =[ and vice versa. Here 
and in the following, the sums over a run over a = ±1. 
Furthermore, we have introduced the abbreviations n™ — 
(l-a)/2 + an a and Eg = E a + [(a + 1)/2]U. The next 
contribution to Eq. (|1L)|) is the renormalization of the dot 
Hamiltonian Ho, 



\Tikf 



Ikaa 



(12) 



As this contribution can be absorbed in the definition of 
Hq, it will be ignored in the following. The last term, 
which describes processes that change the electron num- 
ber on the dot by two, reads 



H2e = 



'kk'. 



Qk - E° 



d s d a 4 k /=cL_ + h.c. (13) 



Later, we shall focus on the sector of the Hilbert space 
where the dot is singly occupied, and hence H 2e can be 
neglected, as well. 

Similarly, we can decompose the phonon-mediated con- 
tributions to the transformed Hamiltonian @ into 

^[L^[L^H T ,H S },H T ] = H$ ex + Hg n + H^ . (14) 

Here, the relevant contribution containing direct and ex- 
change terms reads 

1 
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(15) 



where we have introduced the phonon operator 



Mka a - 2^ — _ pa 
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££k - E£ + frw c 



(16) 



Furthermore, the transformation leads to a renormaliza- 
tion of the spin-phonon coupling 

ffr P en = ~\ E l^k| 2 ^k CT -4^ + h.C. (17) 
£k<r 

Finally, a two-electron tunneling term 
•#2e = 9 E a Tikli/v -X"^ kcra d CT d ff cJ, k/(J cJ k(T + h.c. 



'kk', 



(18) 

is generated, which, for the same reasons as above, will be 
disregarded together with the renormalization term l|17fl . 



IV. CURRENT 

The mean current h{t) across contact I — L, R is given 
by the expectation value of the time-derivative of the to- 
tal number of electrons Ng — X^ko- n £ka- in lead £ multi- 
plied by the electron charge — e, i.e., 



h{t) = -e(N t ) 



: ([H,N e }). 



(19) 



For the evaluation of this expectation value, we switch to 
the interaction picture to obtain 29 - 30 



t-t 



W) - f 



dr<[JVi(t-to),Hi(t 



*o)]>- (20) 



Here, the tilde denotes interaction picture operators with 
respect to H and Hi :— H — H is the coupling part of 
the effective Hamiltonian and the expectation value has 
to be taken at time t — t — to. 

In the long-time limit t 3> to, the mean value of the 
current becomes time-independent and we find that due 
to charge conservation 7l = —Jr. Hence, we can write 
the total current as I = Irl — Ilr, where I in is the con- 
tribution of the electrons flowing from lead i into lead £'. 
In the cotunneling regime, where E a <§; m <C E a + U 
for a =|, I and I = L, R, the dot is, up to exponentially 
small corrections, always occupied by a single electron 
and hence is characterized by the occupation probabili- 
ties p a = (n a ) with P]- + Pi = 1. The currents Ii>i can 
then be written a o 23 i 31 



Iff = e 



(21) 



where We' a >e a is the rate for an electron tunneling from 
lead I into lead £' when the dot was initially in the state a 
and is in the state a' afterwards. By convention, a pro- 
cess with a' = a is called elastic cotunneling, while the 
case with a' — a is referred to as inelastic cotunneling. 22 
The contribution due to the cotunneling without par- 
ticipation of phonons to the rates Wi> a t£ a will be desig- 
nated by Wyl, lcr . In addition, the phonon- assisted co- 
tunneling process contributes with rates W^, (a to the 



4 



total rates Wi> a 'i<j = W^, (cr + W^, ecr . In the next 
two subsection, we discuss these two contributions in de- 
tail. Note that we thereby restrict ourselves to temper- 
atures above the Kondo temperature, where correlations 
between the dot spin and the collective spin of the elec- 
trons in the lead do not play a role. We shall come back 
to a brief discussion of the low-temperature behavior in 
Sec. lWII 



A. Elastic and inelastic cotunneling 

Evaluating the current expression l|2l)[) for the direct 
and exchange terms we arrive at the elastic cotun- 
neling rate 



w: 
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and ne/g — [i^ — \xg. Here, terms of the order 
0( [max(fc B T, \(ig>t + A a > a \) / A^ a > ta ] ) have been ne- 
glected in the curly brackets. 



B. Phonon-assisted elastic and inelastic 
cotunneling 

We now come to the current contributions due to the 
presence of the spin-phonon coupling. Using if| ir cx from 
Eq. I]15p. we obtain after a straightforward but lengthy 
calculation the rates for the phonon-assisted elastic co- 
tunneling 



T -g^ ,deT,,(e') /der,(e)A«( e , e ') 



xN(e-e>)Me) [1 - MO] 



(22) The corresponding energy denominator 



(27) 



with Ai 0) ' cl (e) = {E 9 + U-ef 
the inelastic cotunneling rate is given by 



e) 2 . Similarly, 
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(23) 
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where AS ""^(e) 



-n 2 



{E ff + U-£) l -{E a -e) 1 \\ In 
the inelastic case, the energy A Sa is deposited in or ex- 
tracted from the dot. Correspondingly, the initial en- 
ergy e of the tunneling electron coming from lead I differs 
by the same amount from the final energy of the elec- 
tron in lead £'. Inelastic cotunneling is, thus, strongly 
suppressed for temperatures and external voltages much 
smaller than A Sa . 

To a very good approximation, the energy dependence 
of the coupling rates I^(e) can be ignored in the en- 
ergy range contributing to the integrals in the rate ex- 
pressions <|22ll and l|23|) . Furthermore, near one of the 
resonances, e.g., for E a < m -C E„ + U, we can ex- 
pand the dominant energy denominator of A^ (e) and 
A^' inel (e), respectively, and carry out the energy inte- 
gration^ This yields the approximate rates 



W {0) 



fiTiJV 1 



2tt 



A 2 



1 



4A 2 n 



47r 2 (fc B T) 2 + (n eu + A al(j y \Q{-mu - A a , a ) 



where we have introduced the function 

e 



6(e) 



1 - exp(-e/fc B T) 



as well as the energies 



E„ 



E„ 



(24) 



(25) 



(26) 



ai: , (^')=e 



(E a + U- e) (E s + U- e') 



+ 



(E„ - e) {E s - e') 



(28) 



now depends on both the energy before and after the 
phonon emission and/or absorption. Note that as ex- 
plained above, we use the term "elastic" only to indicate 
that the dot state before and after the cotunneling pro- 
cess is identical. From the point of view of the electron 
system, the phonon-mediated cotunneling process is, of 
course, no longer elastic. Rather, the electrons either 
emit the energy e — e' due to the stimulated and spon- 
taneous emission of phonons or absorb the energy e' — e 
from the phonon system. Both processes are captured by 
the function 



N(Ae) = D(Ae/h) [n B (Ae) + 1] + D{-Ae/K) n B (-Ae) . 

Similarly, we obtain the "inelastic" phonon-mediated 
cotunneling rates 



W vL = ^ /de'l>(e') |der,( e )A«(e,e') 

xN{e-e' -A 9 „)f t {e) [1 - MO 

with 

1 



(30) 



(E a + U-e)(E (7 + U-e') 



(E a - e){E a - e>) 



(31) 



For the further evaluation of the rate formulas we again 
assume energy- independent dot-lead couplings Tp. Since 
the energy denominators in Eqs. (|28|l and l|31|l are al- 
ready of fourth order in A CTAI , we can neglect their energy 
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dependence around e = e' = p, to get a result valid to the 
same order as in Eq. \2ty . Carrying out one of the en- 
ergy integrals, we obtain the approximate result for both 
rates 



(!) KTiTt, 
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8tt 2 V A 2 A 2 A 2 A 2 



x /dAeAr(Ae)6(- m -A CT v-Ae) (32) 

with the energy differences A a £ = E a — /i£. We now 
consider a spectral density of the spin-phonon coupling 
of the form 



D{uj > 0) = jhujl- s uj s 



(33) 



with the dimensionless spin-phonon coupling parame- 
ter 7, the cut-off frequency lu c and an exponent s, which 
depends on the nature of the coupling. In order for our 
master-equation approach to be valid, we have to re- 
strict ourselves to exponents s > 1. The exponent s = 1 
corresponds to the so-called Ohmic case, and below, we 
will find s = 3 or s = 5 for a spin-orbit mediated cou- 
pling to phonons. Inserting the spectral density i|33|) into 
Eq. 132|) , we can write the rates as 



rxr(l) lb \l-s^I> / 1 



i 



8tt 2 V A 2 A 2 A 2 A 2 

x (fc B r) 1+s J, (- (m + a^ ct ) /k B T) e (- m - a ctV ) 



(34) 



with 



l-e~ a /"°°, „ N |lnd s lnx + a , . 

J s (a) = / dx sgn(lna: J — . (35) 

a j x — 1 x — e 

In the interesting cases of s = 1, 3 and 5, we find for the 
integral 



J 1 (a) = ±(4 n * + a 2 ) 

J 3 (a) = i_(4^ 2 +a 2 ) (8^ 2 + 3a 2 ) , 



(36) 
(37) 



J 5 (a) = — (2tt 2 + a 2 ) (4tt 2 + a 2 ) (8vr 2 + a 2 ) . (38) 

In the general case, one can readily show that J s (a) is an 
even function, i.e., J s (a) — J s (—a) for all s and a. 



V. MASTER EQUATION 

For the evaluation of the current expressions given in 
the previous section, we still need the occupation prob- 
abilities p a . Their dynamics is governed by the master 
equation^ 



W- ot 

rV (7(7 



W?i P . Here, Wg* 



with the total rates W aa 
J2w Wi'gia is the spin-flip rate due to inelastic cotun- 
neling, which, in contrast to the current formula 121fl . 
also contains processes involving a single lead only. Fur- 
thermore, spin-flip processes due to the spin-phonon cou- 
pling (jSJ) have to be taken into account in the master 
equation l|39|l . These processes lead to the additional 
rate 



W^ = -N{K S ) 



(40) 



In the stationary limit, the solution of the master equa- 
tion is given by p a = [1 + Wg a /W a g\~ , where we have 
used the normalization condition P] + Pi = 1. Note that 
since elastic cotunneling leaves the dot state invariant, 
besides the spin-flip rate W^ lp only the rates for the in- 
elastic cotunneling processes appear in the master equa- 
tion. 

Figure [21 shows the ratio P]/pi of the populations as 
a function of the bias voltage V in the presence of an 
Ohmic environment for different coupling strengths 7. 
Near thermal equilibrium, i.e., for voltages far below the 
onset of inelastic cotunneling, i.e., eV <C |Ajj|, k B T, the 
total rates and thus the populations fulfil the detailed 
balance relation p s /p a — W Sa /W aa = exp(—A Sr7 /k B T). 
For higher bias voltages, inelastic cotunneling leads to a 
population of the excited dot state, i.e., a heating of the 
dot. This heating effect becomes less pronounced in the 
presence of the spin-phonon coupling, where the relax- 
ation due to the rate (|40|l equilibrates the system. For 
a quantitative estimate of the heating effect, we consider 
the deviation h of the population ratio from its equilib- 
rium value: 



h = 



Pi 



(41) 



Using the detailed balance relation for the rates Wg 
we then find 



flip 

(7(7 5 



Pa = -W Sa p a + W aa p a 



(39) 



h _ |exp(A iT /fc B r)-ITff/^f| 

i + wffW 

Thus, even for a non-zero difference in the numerator, 
the heating effect becomes suppressed for W^ p ^> Wm?* . 
In Fig. [2J we indicate the onset voltages of the heating 
regime, defined by the relation W^ lp = WPf* by vertical 
dashed lines. 

Finally, we remark that heating effects are suppressed 
when the dot-lead coupling is asymmetric^ say Tl 3> 
Tr. In this case the excitation of the dot, which is domi- 
nated by processes involving the tunneling from one lead 
to the other and which, hence, is proportional to the 
product TiTr , is negligible compared to the relaxing ef- 
fects resulting from the coupling to the left lead, which 
are proportional to r£. Such a situation may occur, e.g., 
when measuring the current through an atom or molecule 
on a surface which is contacted on one side by an STM 
tipiS (cf. also discussion at the end of Sec. IVIl| l. 
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FIG. 2: Ratio /py of the populations as a function of the 
bias voltage V in the presence of an Ohmic bath with different 
coupling strengths 7. The parameters are Ajf = O.lEo, T = 
O.OlEo and KT L = ftF R = 0.02E . Here, E = \Ei - n\ is the 
energy difference between the energy E± — E-j — Ey of the 
degenerate dot level and the Fermi energy fi = fit, = ftR in 
the leads in the absence of both the magnetic field and the 
bias voltage. The dashed vertical lines indicate the onset of 
the heating regime defined by W+? p 



VI. DIFFERENTIAL CONDUCTANCE 

Having all the information for the evaluation of the 
current l|21|l at hand, we can now turn to the experimen- 
tally relevant quantity, the differential conductance g = 
dl/dV. As a function of the bias voltage, this quan- 
tity shows a characteristic step at the onset of inelas- 
tic cotunneling around |A|||/e, when another transport 
channel becomes accessible. It is known that the width 
of this step is proportional to the temperature T, while 



the broadening of the dot states due to their coupling to 
the leads does not play a role. Compared to a measure- 
ment in the sequential tunneling regime, an experiment 
can thus achieve a much higher precision for the deter- 
mination of the splitting A^ (see Refs. Hoi and fl2|) . An 
open question, however, is how spin-flip processes taking 
place during the cotunneling process affect such a mea- 
surement. 

Let us first consider the zero-temperature case with- 
out heating where the dot is always in its ground state 
which we assume to be the spin-up state. Evaluating 
the current using the approximate cotunneling rates l|24|) 
and i|34|) , we then find for the contribution due to elastic 
cotunneling 
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[4TT 2 (k B T) 2 + 3(eV) 2 



47T 



' A 2 
LTRT 

A2 A2 



A 2 

A2 A2 



x (k B T) s [k B TJ s (eV/k B T) + eVJ' s (eV/AfeT)] } , (43) 



where J' s (a) = dJ s (a)/da. In order to obtain a compact 
expression for the inelastic cotunneling current, we fur- 
thermore assume positive voltages V > and not too 
high temperatures k B T <C Ajj. Then only the tunneling 
from the left to the right contact is relevant, i.e., the total 
current / is to a very good approximation given by Trx- 
Furthermore, away from the step around V ~ Ajf/e, 
we can use that up to exponentially small corrections 
6(e) rs for e < k B T and 6(e) » e for e > k B T. Thus, 
for eV <C A||, the inelastic cotunneling current vanishes 
and for voltages eV ^S> A^j, it contributes with 



fed = + W 2 (k B Tf + 3(eV- A iT ) 2 " 

+ ^^f 1 ^ ( + 3^r) (W [k B TJ s ((eV - A iT ) /kBT) + (eV - A iT ) J' s {(eV A iT ) /k B T) } ) 

(44) 



to the differential conductance. Particularly interesting is 
the step in between, which is centered around V = A^ / e. 
In order to obtain an estimate for its width AV s t C p, we 
determine the onset V on of inelastic cotunneling by us- 
ing a linear approximation of the differential conduc- 
tance around the center of the step, g(V) sw g{A^/e) + 
g'(A^/e) (V — Ajj/e) and intersecting with the purely 
elastic cotunneling curve lf4l5|l (see Fig. OJ: g c \{V n) — 
5(A;t/e) +g'(A iT /e) (V on - A iT /e). Approximately, we 



can determine the solution of this non-linear equation by 
employing a linear approximation for the elastic cotun- 
neling curve around the center of the step. Doing so, we 
obtain V on = A iT /e-5 in oi(A iT /e)/g( nel (A 1T /e). Defining 
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FIG. 3: Sketch of the procedure employed for the determina- 
tion of the onset Von and of the width AV a t ep of the inelastic- 
cotunneling step. The dashed line shows the extrapolation of 
the differential conductance curve below the step. The dash- 
dotted line corresponds to the linear approximation of the 
step shape. 



the width as twice the difference V on — Aj^/e gives 



AV stC p = 



6k B T 



1 - 3 



fc R T 



2 r 



1 + 7 



1 

6tt 



A 4 

A 2 A 2 



A 4 



A 2 A 2 



x (3J' S '(0) + 2J' S (0)) 



hhJr . 



(45) 



As mentioned above, the width is essentially propor- 
tional to the temperature. We also find that the spin- 
flip processes only modify the higher-order corrections in 
/cbT/At4R|. Here, the spin-flip coupling contribution is 
proportional to 7 and to a positive prefactor of the order 
of one and has a temperature dependence which depends 
on the ratio of thermal and cutoff energy. In fact, we 
find that with increasing coupling constant 7 the width 
of the inelastic cotunneling step becomes reduced. 

We now turn to the discussion of the heating effects 
by taking into account the stationary probability distri- 
bution according to the master equation fl35fr. Figure 0] 
shows the differential conductance as a function of the 
voltage in the presence of an Ohmic spin-flip coupling of 
strength 7. In the upper panel, the case of a symmetric 
dot-lead coupling is depicted. As discussed above, this is 
the situation where heating effects are most relevant, and 
for 7 = we indeed observe a noticeable deviation from 
the step-like behavior for voltages V > A^/e. With in- 
creasing spin-flip coupling strength 7, however, the dot 
is driven to equilibrium and the overshooting of the dif- 
ferential conductance curve disappears. As mentioned at 
the end of Sec. E] heating effects are suppressed when 
the dot-lead coupling is strongly asymmetric. Such a sit- 
uation is shown in panel (b) of Fig. In this case, the 
influence of the spin-phonon coupling is much weaker. 



For larger voltages, the effect of phonon-assisted inelas- 
tic cotunneling becomes visible, as can be inferred from 
the inset of panel (b) . 
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FIG. 4: Cotunneling conductance g as a function of the bias 
voltage V in the presence of an Ohmic bath with different 
coupling strengths 7. The parameters are A|t = 0.1-Bo, T = 
O.Ol^o and (a) ftF L = hV R = 0.02.Eo, (b) fiT L = 0.2E and 
HTl = 0.002i?o with Eo as defined in the caption of Fig. |5] 
In the inset of panel (b), the relative difference between the 
differential conductance g in the presence of the spin-phonon 
coupling and the differential conductance go without spin- 
phonon coupling is depicted. 

In Fig. [3 the step-width defined by the onset of in- 
elastic cotunneling is shown as a function of the spin-flip 
coupling strength 7 for different ratios T^/Tr. For equal 
coupling to both leads, the main effect again comes from 
the suppression of the heating with increasing 7. The in- 
set of Fig.0depicts the 7-dependence for larger 7. As de- 
scribed by the analytical expression (|45fl the width drops 
linearly. Due to the finite temperature, there is a slight 
offset compared to the full solution, however. 



VII. COMPARISON WITH EXPERIMENTS 

We now apply the general results to the specific ex- 
perimental situation of Ref. [12| and analyze the rele- 
vance of spin-phonon coupling effects in this case. In 
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0.595 




2.5 x 10 



7.5 x 10" 



1 x 10~ 3 



FIG. 5: Step width AV stcp = 2 ft ' ncl (A iT /e)/g inc i(A u /e) as 
a function of the coupling strength 7 to an Ohmic bath for 
different ratios of the tunneling rates Tl.r, where the product 
h^TiTn = (0.02_Eo) 2 is kept constant. The other parameters 
are like in Fig. 2] The horizontal dashed line shows the width 
according to Eq. ^451 for 7 = 0. In the inset, the behavior 
for larger 7 is shown together with the step width without 
heating according to Eq. 1451 (dashed). 



order to do so, we model the 2D quantum dot by an 
isotropic parabolic lateral confinement characterized by 
the frequency loq in the absence of a magnetic field. In a 
magnetic field B z perpendicular to the two-dimensional 
electron gas (2DEG), the relevant frequencies then are 
u>x,% = fi T w c /2 with the cyclotron frequency to c = 
\e\B z /m*c, where m* is the effective electron mass and 
ft = y/ojQ + l^c /4 corresponding to the effective lateral 
confinement length I = yjhm* /ft (for more details, we 
refer the reader to R,e f.ll8fi . The parameters for the GaAs 
quantum dot of Ref. Il2l are g = —0.16, m* — 0.067m c , 
B z = 11T and we assume Hluq — l.lmeV. 

We consider the influence of both Rashba^ and Dres- 
selhaus 34 spin-orbit coupling, which are characterized by 
spin-orbit lengths Ar = Ad = 8/im (see Appendix A). 
For a piezoelectric coupling to longitudinal and transver- 
sal phonon modes in the limit u> <C c a min(f, d), where 
d w 5nm is the vertical confinement length, the spectral 
density (|33[1 of the spin-phonon coupling is cubic, i.e., 
s = spe = 3. The prefactors for the longitudinal and 
transverse contributions are given by 



7pe,l 



7pe,t 



Qb + g| ( eh u 
70nhgc? 



Ql + Q 



1057r?i gc^ 



(46) 
(47) 



with the speed of sound cl = 4.73 x 10 5 cm/s and 
ct = 3.35 x 10 5 cm/s of the longitudinal and transver- 
sal acoustic phonon modes, respectively. The density 
of GaAs is g — 5.319 g/cm 3 and the electron-phonon 
coupling parameters are eh%4 = 1.2 x 10 7 eV/cm and 
k = 13.2. The spin-orbit coupling strength enters via 



the dimensionless coupling constantsiZii& 



Qr,d — 



I 



R,D 



Hu>i 



hu>9 



Huji =F Aj_j hiL>2 ± A 



11 



(48) 



The cut-off frequencies for the piezo-electric coupling are 
given by w Ci pe,l = cl/1 and w c ,pe,t = ct/1. Note that 
the leading contribution to Qr,d is linear in B z , and 
consequently the rate l|40(l is proportional to B\ for low 
magnetic fields B z . 

In the case of a deformation potential coupling to 
acoustic phonons, the spectral density is proportional 
to w b , i.e., s = sda = 5, and the prefactor is given by 



7DA 



^ (Q 



Ql) 



2Anh g I 



2 „3 



(49) 



with S = 6.7 eV. The cut-off is at the frequency w Ci da = 
cl//. Taking into account the magnetic- field-dependence 
of the prefactor 7da, the leading contribution to the 
rate (f4"0")l is proportional to B 7 Z at low magnetic fields. 
Thus, the spin-phonon coupling due to the deformation 
potential electron-phonon interaction is only relevant for 
rather high magnetic fields. Furthermore, we take a 
temperature T = 15 mK and dot-lead coupling rates 
hT L = HT U = 17.5 ^eV from Ref. O Finally, we as- 
sume a value \E\ — = 0.4meV, where E\ and fj, are 
the energy of the degenerate dot-level and the Fermi en- 
ergy, respectively, in the absence of the bias voltage and 
the magnetic field. An external bias is assumed to drop 
symmetrically across both contacts. 

Figure|fj]depicts the differential conductance as a func- 
tion of the bias voltage for the given set of parameters. 
Comparing the result of the calculation including the 
spin-orbit induced spin-phonon coupling (solid line) with 
the one without this coupling (crosses), we can conclude 
that these effects play a negligible role in the experiment 
by Kogan et al. However, as shown in the inset of Fig. El 
by reducing the dot-lead coupling by a bit more than 
one order of magnitude, one reaches a regime where the 
equilibration of the dot due to the spin-phonon coupling 
becomes relevant. As discussed in Sec. the transition 
to this regime is determined by the relation W^! p = 
which for eV > Aj,j yields the tunneling rate 



^cot 

' IT ' 



on, he at 



= IA 



2ttW. 



URTl 



flip 
11 



(50) 



This relation can be employed for the determination of 
the spin-flip rate W^ p by means of a cotunneling trans- 
port measurement. In order to do so, one has to de- 
termine the critical tunneling rate r onj hcat marking the 
transition between the regime of the dot being in equilib- 
rium and the one where heating effects are relevant. As 
an indicator for this transition one can measure the ratio 
of the differential conductance above and below the in- 
elastic cotunneling step. Determining this quantity as a 
function of the dot-lead coupling strength yields r onj h ea t 
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FIG. 6: Differential conductance g as a function of the bias 
voltage V for the parameters given in the main text. The 
crosses show the result in the absence of the spin-phonon cou- 
pling. In the inset, the ratio of the differential conductance 
gi at Vi = 1.5A||/e (see arrow in the main panel) and the 
differential conductance go at zero voltage for different values 
of the dot-lead coupling rates V = Tl = Tr. The dashed 
vertical line designates the value given by Eq. 1501 and the 
arrow the dot-lead coupling rate in the experiment—. 



(cf. inset of Fig. [fJJ) and via Eq. (|50l) the spin-flip rate 



flip 

1! 



^1 

Let us finally discuss whether such a measurement 
scheme would also apply to transport experiments 
through single atoms or molecules on a surface, e.g., the 
ones reported in Ref. There, the tunnel-coupling to 
the two leads is strongly asymmetric, one contact being 
via the insulating surface and the other via an STM tip. 
As has been discussed at the end of Sec. |Vj heating ef- 
fects are suppressed in this case. Indeed, the form of 
the differential conductance-voltage characteristics seen 
in Ref. llfj can — above the Kondo temperature — be very 
well described by an equilibrium probability distribution 
of the spin state of the atom. However, for a not too small 
tunnel-coupling across the insulating layer, it should be 
possible to enter the heating regime by adjusting the 
atom-tip coupling strength appropriately. If this cou- 
pling strength is not too far away from the critical value 
given by Eq. J5UJ), it should be possible to measure the 
intrinsic spin-relaxation rate due to phonons by using a 
similar scheme as the one proposed above, although only 
the atom-STM coupling strength can be controlled easily. 
On the other hand, it is in principle feasible to control 
the molecule-surface coupling strength by modifying the 
chemical binding of the molecule to the surface, as was 
shown in Ref. Illl so the proposed scheme could be ap- 
plied to such a situation, as well. 



VIII. SCALING BEHAVIOR 



temperature Tk- For lower temperatures, correlations 
between between the dot spin and a collective spin degree 
of freedom of the electrons in the leads become impor- 
tant and a perturbative treatment of the problem breaks 
down. This leads to the so-called Kondo effect, which 
has been analyzed in great detail in the literature for 
the Anderson models. We now briefly discuss the ef- 
fect of additional spin-flip processes due to the term Hg 
using the poor man's scaling approach put forward by 
Anderson^. In order to keep the discussion as simple as 
possible, we consider a simplified version of our effective 
Hamiltonian @, which ignores the lead-, momentum- 
and spin-dependence of the coupling constants as well as 
the direct interaction between the lead states: 

H = H - JS d • s + J s (XS% + YS v d ) 

-J 1 (XSl + YSl)n(0) + J 2 (Xs x + Y s y) { ' 

Here, we have introduced the spin operators 

s = - ^2 4ka T °°' c e'k'a' and ( 52 ) 

«'kkW 

Sd = - <4 Too' da' , (53) 

where r is the vector consisting of the three Pauli matri- 
ces. The density of lead electrons at the position of the 
dot is given by 



(54) 



ee'kk'a 



The phonon operators are defined as X = ^ q M^ x (a q + 

a f _ q ) and Y = £ q M q>!/ (a q + al q ). 

During the scaling procedure, contributions due to ex- 
citations of electrons and holes far away from the Fermi 
energy = 0, at a cutoff energy ±D, are successively 
integrated out. If the resulting effective interaction is 
again of the form (|51|l . we can absorb it in the coupling 
constants. Note that in general, the coupling constants 
thereby acquire a momentum dependence. However, as 
long as the cutoff energy is still large enough, this de- 
pendence can be ignored for the relevant states near the 
Fermi energy. Doing so, we obtain the usual scaling equa- 
tions for the exchange constant J: 



dJ_ 
dD 



-2v 



J2 
D 



(55) 



which, in particular, is not influenced by the other terms 
in the Hamiltonian (|5f I) . The renormalization of the 
electron-phonon coupling constant is given by 



dJs TT le ^ D 1 

^- vJJ2 \ FTw 



(4 In 2) v 2 JJ 2 (56) 



So far, we have restricted the discussion to the 
case of temperatures well above the so-called Kondo- 



while the coupling constants J\ and J 2 remain invari- 
ant. Here, v is the density of states in the leads. In 
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the anti-ferromagnetic case, where J > scales to oo as 
D — ► 0, the same holds true for the spin-phonon cou- 
pling constant Jg. We conclude that far away from the 
active region around the Fermi energy, neither the pure 
spin-flip term proportional to Jg nor the phonon-assisted 
tunneling terms, i.e., the last two terms in the Hamilto- 
nian (|51|l , modify the scaling behavior of the Kondo-type 
coupling of the dot spin to the leads. 



IX. CONCLUSIONS 

We have theoretically studied the influence of a spin- 
phonon coupling on the cotunneling through a nanoscale 
system. By means of a Schrieffer- Wolff transformation, 
we have derived an explicit Hamiltonian in which the rel- 
evant cotunneling processes appear to lowest order. This 
allowed us to evaluate the rates for elastic and inelastic 
cotunneling in the presence of spin-flip processes. We 
found that the width of the inelastic cotunneling step is 
only weakly influenced by the spin- flip processes. More 
important are the relaxation effects of the spin-phonon 
coupling which counteract the heating of the quantum 
dot due to inelastic cotunneling. Considering a realistic 
model for a recent experiment 12 , we propose a new way of 
determining the spin relaxation rate due to spin-phonon 
coupling in a cotunneling experiment. 
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APPENDIX A: DERIVATION OF SPIN-PHONON 
COUPLING HAMILTONIAN 

In this appendix, we derive the spin-phonon cou- 
pling (JjJJ starting from a Hamiltonian describing an elec- 
tron in a two-dimensional electron gas with an external 
magnetic field B z in z-dircction 



Ho 



p2 1 



2m 



(Al) 



where P = p + (|e|/c)A(r) is the kinetic momen- 
tum, U(r) is the confinement potential, and A(r) = 
(B z /2)(—y, x, 0) is the vector potential. For x, y, z point- 
ing along the main crystallographic axes of the GaAs 
crystal, and the electron gas lying in the [001] plane, the 
spin-orbit interaction assumes the form 



with the Dresselhaus^ and RashbaS contributions 



Hb = (3n{-<J x P x + (JyPy) , #r = a R (a x P y - a y P x ) . 

The strength of the spin-orbit coupling is conveniently 
measured in terms of the spin-orbit lengths Ad = 
h/m* (3-q and Ar = fi/m*aR. Finally, the orbital degrees 
of freedom of the electron are coupled to the phonon 
system via the electron-phonon interaction i/ _ p h = 

Eq-^qWK + a -q)> where a branch index has been 
suppressed. 

Applying to the Hamiltonian H = ZZo + i?so+-He-ph a 
Schrieffer- Wolff transformation which eliminates Hgo to 
lowest order and projecting on the orbital ground state, 
we obtain to lowest order the spin-phonon coupling 



H s = (0|[5,fle-ph]|0) 



(A4) 



where the generator S of the canonical transformation 
has to fulfil [S,H ] = -H SO - From Eqs. (SlHIsSl, we 
can infer that the generator can be written in the form 
S = i(A x a x + AyO~y) with A x and A y being Hermitian 
operators which act on the orbital degrees of freedom 
only. Hence, we obtain 

H s = i(0|[A X)J ff e _ ph ]|0) ( r*+i(0p 3/)J ff e _ ph ]|0) a y . 

(A5) 

Letting M q ,j = (0|[f2j,M q (r)]|0) for i = x,y, we thus 
arrive at the coupling Hamiltonian JSJ. 

For a quantum dot with parabolic lateral confinement 
(see Sec. Will , it is useful to choose coordinates rotated 
by 45° degrees around the z-axisii, i.e., we let x — > (x + 
y)/v2 and y — > (y — x)j\/2. In this new coordinate 
frame, the generator S of the canonical transformation is 
defined in terms of 
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LU 2 



huJ2 — Ajf 

— m*Clx — p y 

>2 T 7 

nw2 + A|| 
m*Qx + p y 



hw2 — A 



LT 



(A6) 



(A7) 



H S o = H B + H r 



(A2) 



Note that this transformation diverges for magnetic field 
strengths with hwi.2 ± Ajf = 0. A detailed analysis of 
the relaxation rates in this case can be found in Ref. 0. 
On the other hand, for \A^\ <C ^1,2, the leading con- 
tribution to the relaxation rates comes from the terms 
linear in Aj_| cx B z (see Ref. ^ 
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